arXiv:hep-th/9502011vl 2 Feb 1995 


USC-95-01 

hep-th/9502011 


Boundary bound states and boundary bootstrap 
in the sine-Gordon model 
with Dirichlet boundary conditions 


Sergei Skorik^ and Hubert SaleuD* 

'Department of Physics, University of Southern California 
Los Angeles CA 90089-0484 

*Department of Physics and Department of Mathematics 
University of Southern California 
Los Angeles CA 90089 

We present a complete study of boundary bound states and related boundary S- 
matrices for the sine-Gordon model with Dirichlet boundary conditions. Our approach 
is based partly on the bootstrap procedure, and partly on the explicit solution of the 
inhomogeneous XXZ model with boundary magnetic field and of the boundary Thirring 
model. We identify boundary bound states with new “boundary strings” in the Bethe 
ansatz. The boundary energy is also computed. 


* Packard Fellow 


February, 1995 


1. Introduction. 


The sine-Gordon model with a boundary interaction preserving integrability (which 
we shall call the boundary sine-Gordon model) is of theoretical as well as practical interest. 
In particular, it exhibits relations with the theory of Jack symmetric functions [jT| and has 
appplications to dissipative quantum mechanics J2) and impurity problems in ID strongly 
correlated electron gas |§ . 

In the seminal work ||| it appeared clearly that this problem presents an extremely 
rich structure of boundary bound states, which was partly explored in ||. Our first 
purpose here is to study this structure further in the particular case of Dirichlet boundary 
conditions, that is the model 


£sg = - 


(dtV? 


m r 


(d x v ) 2 + ^cos fa 


dx 


( 1 . 1 ) 


with a fixed value of the field at the boundary: (p(x — 0, t) — ipo. 

Also, the consideration of boundary problems poses interesting challenges from the 
point of view of lattice models, here lattice regularizations of ( |1.1|) . In |6] and also in 
[7j] it was shown in particular how to derive the S-matrices of |4|| from the Betlie ansatz. 
Our second purpose is to complete these studies by investigating which new types of 
strings correspond to boundary bound states, and by deriving as well the set of S-matrices 
necessary to close the bootstrap. Observe that lattice regularizations are useful to define 
what one means by putting a bound state at the boundary. Indeed, some bound states 
have no straightforward interpretation, and although they are easy to study formally using 
the Yang Baxter equation and the bootstrap, their meaning in the field theory is unclear. 

In section 2 we consider the bootstrap problem directly in the continuum theory. We 
identify boundary bound states and we compute the related boundary S matrices. In 
section 3 we write the Betlie ansatz equations for the inhomogeneous six-vertex model 
with boundary magnetic field, which is believed [6| to be a regularization of o- We 
show that these equations are also the bare equations for the Thirring model with 17(1)- 
preserving boundary interaction, which is the fermionized version of ( |1.1| ). In section 4 
we discuss in details new solutions (“boundary strings”) to the Bethe ansatz equations 
made possible by the appearance of boundary terms. In section 5 we study the physical 
properties of the model, in particular the masses and S-matrices corresponding to these 
boundary strings, and we partially complete the identification with the bootstrap results of 
section 2. Several remarks, in particular formula for the boundary energy of the boundary 
sine-Gordon model, are collected in the conclusion. 
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2. Boundary bootstrap results. 

2.1. Solving the boundary bootstrap equations. 

The S-matrices for the scattering of a soliton (P + ) and an anti-soliton (P _ ) on the 
ground state |0)s of the sine-Gordon model with Dirichlet boundary conditions ( |1 . 1|) were 
obtained in ||]: 

P ± {9) =cos(£±\u)R 0 (u)R 1 (u,£), (2.1) 


where 9 — iu is the rapidity, £ = and A = jp — 1. The explicit form of Ro,Ri 

is rather cumbersome and can be found in jp. Since the theory is invariant under the 
simultaneous transformations £ —> —£, and soliton—>anti-soliton , we choose hereafter £ to 
be a generic number in the interval 0 < £ < Atx 2 / f3 2 (about the value of the upper bound 
see the discussion in ||). 

The function Rq contains poles in the physical strip 0 < Im 9 < tt/ 2 located at 
u = n = 1, 2,... < A. These poles arise because of the corresponding breather pole in 
the soliton-antisoliton bulk scattering, and should not be interpreted as boundary bound 
states 0. 

When £ > n/2, the function P + (9) has additional poles in the physical strip , located 
at u — v n with : 


0 < v n 


£ 

A 


2n + 1 
2A 


7T < 


7T 

2 ’ 


( 2 . 2 ) 


(n = 0, 1, 2,...) corresponding to a first set of boundary bound states which we denote by 
| f3 n ), with masses 


m 


n 


m cos v n = m cos 



2n + 1 A 
~~2X~ 7T ) ’ 


(2.3) 


where m is the soliton mass. These bound states are easy to interpret For 0 < 

ifo < the ground state of the theory is characterized by the asymptotic behaviour tp — > 0 
as x —> oc. but other states, whose energy differs from the ground state by a boundary 
term only, can be obtained with <p —> a multiple of ^ as x —> oo. Since the f3 n appear 
as bound states for soliton scattering, they all have the same topological charge as the 
soliton, which we take equal to unity by convention, so they are all associated with the 
same classical solution, a soliton sitting next to the boundary and performing a motion 
periodic in time (”breathing”), with — 0) = ipo and <p —> j as i -> oo 0. 

To deduce the scattering matrices on the boundary bound states we use the ’’boundary 
bootstrap equations” as given in |ij. We assume that these S-matrices are diagonal, which 
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is true if all the boundary bound states have different energies. In this case the bootstrap 
equations read: 

rtjs W = E + <a)S£(« - <,)■ (2.4) 

c,d 

These equations allow us to find the scattering matrix of any particle b on the boundary 
bound state /? provided that the latter appears as a virtual state in the scattering of the 
particle a on the boundary state a. The masses of the corresponding boundary states are 
related through 

mp = m a + m a cosuf a , (2.5) 


where iv^ a denotes the position of the pole, corresponding to the bound state (3. 

Let /3 n stand for the n-tli boundary bound state corresponding to the pole v n in P + 
(|2.2|). Then (|2.4| ) gives: 


Pp n (9) = P + (9)a(9 - iv n )a(9 + iv n ), 


( 2 . 6 ) 


Pp n {0) = b{6-iv n )b{6+ iv n )P (9) + c{9 - iv n )c(9 + iv n )P + (9), (2.7) 

where the well known bulk S-matrix elements a(9) — 5++ = SZZ (kink-kink scattering), 
b(9) — S+Z = <ST+ (kink-anti-kink transmission) and c(9) — ST_ = S_+ (kink-anti-kink 
reflection) can be found in [[ji¬ 
lt is easy to check that the matrix elements ( j2.6| )-( |2T| ) satisfy general requirements 
for the boundary S-matrices, such as boundary unitarity and boundary crossing-symmetry 
conditions [fij, e.g. 

PpS y ~e)= b(26)P+(- + 6) + c(20)Pp n (- + 9), 

= 1 - 

Finally we obtain from (|2.6|)-(|T7|) by direct calculation: 



cos(^ - Att - iXO) . . 
cos(£ — A7T + i\6) @ n 


( 2 . 8 ) 


Hence the boundary Yang Baxter equation is satisfied since the ratio of the above two 
amplitudes has a form similar to ( |2 . 1| ) with £ —»> £ — \ir, £ being a free parameter. 

The analytic structure of P^ ( 9 ) is as follows. The function P\ ( 9 ) has simple poles 
in the physical strip located at u — j + 2 ^ 1 7r , N — 0,1,2..., and at u — v n . It has 
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double poles at u — iv n + ik — 1,2, ...n. The function Fy ( 6 ) possesses in the physical 
strip the same singularities as Pp n (0) plus the set of simple poles at u — iw^ with 


WN = 7T - - 


£ 2 AT — 1 , , 1 ^ ^ A + 1 £ 

A+ 2~x >iV> — -x 


(2.9) 


Interpreting these poles in terms of boundary bound states requires some care. First, 
due to the relation (|2.4|), one sees that if (3 appears as a boundary bound state for scattering 
of a on a, then the poles of the amplitude for scattering of b on a are also in general poles 
of the amplitude for scattering of b on f3. It seems unlikely that these poles correspond to 
new bound states, although in our case they would have a natural physical meaning, for 
example one could try to associate them with classical solutions where ip —> ^ as x —> oo. 
Indeed there are strong constraints coming from statistics that we should not forget. For 
instance at the free fermion point /3 1 2 — 47r, there is a bound state (3\. but although /A has 
again a pole at /5i, the state of mass 2 mp 1 is not allowed from Pauli exclusion principle, as 
can easily be checked on the direct solution of the model (see below section 3.3). Therefore 
we take the point of view that the poles already present in the scattering on an “empty 
boundary” are “redundant”. The only poles we interpret as related to new boundary 
bound states are ( |2.9|) (the additional poles in P~ are related to them by crossing). We 
denote these boundary bound states |<5-n, at) , and their masses, according to ( |2.-"{ ) and ( |2.2|) , 
are given by 

, , a (£ 2n+l ^ (i , 2N-1 

Tfl n ^N — m(COSV n + COS W]y) — 171 COS I —-— 7T 1 — m COS I — H-— 7T 

y A J y A 2 A 

b . (£ N-n-l \ 

= m N+n sm + 2A «) , 

where m b p = 2m sin (|^) is the mass of the p-th breather, p — 1, 2,... < A. 

To understand the physical meaning of these new boundary bound states it is helpful 
to consider the semi-classical limit A —> oo of the sine-Gordon model. As discussed above, 
the boundary bound state f3 n , corresponding to ( |2.2j ). are associated to solutions where 
a soliton is sitting next to the boundary and “breathing”. An incoming anti-soliton can 
couple to this soliton, and together they form a breather sitting next to the boundary and 
performing again some (rather complicated) motion periodic in timcS. The quantization 

1 To compute this solution explicitely requires using a bulk five-soliton configuration Q, an 

expression which is very cumbersome. 
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of this solution shoud lead to |<5 nj jv)- The topological charge of the states |is equal 
to 0 in our units, or, equivalently, to the charge of a free breather in the theory ( |1.1|) . 

One can in principle continue to solve the bootstrap equations ( |2.4| ) recursively. For 
example, for the scattering of solitons or antisolitons on the boundary bound states |h n ,.v) 
( |2.9| ) one obtains the following S-matrices: 


P 5 n , N ( e ) = Pp n (0)a(0-iw N )a(0+ iw N ), 


( 2 . 11 ) 


p+ = cos(£ - i\6) p _ 
5 ™’ n cos(f+ i\9) Sn ’ N ' 


( 2 . 12 ) 


Pf~ has only one simple pole in the physical strip at u — Wn, while P$ x has also simple 
poles at u — Vk, k — n + l,n + 2,..., [| — |j. According to the discussion below (pT9| ) , 
we do not consider these poles as associated with new boundary bound states. Therefore, 
the boundary bootstrap is closed for solitons and antisolitons in the sense that further 
recursion will not generate new boundary bound states. 

So far we have obtained two sets of boundary bound states ( | 2 . 3 | ) and ( | 2 . 10 | ) by con¬ 
sidering all the poles in the physical strip of amplitudes for scattering a soliton (resp. 
anti-soliton) on a boundary with or without a boundary bound state. Of course we should 
also consider the scattering of breathers off the boundary. The scattering of breathers off 
an “empty” boundary was studied in ||, and we refer the reader to this work for the ex¬ 
plicit boundary S-matrices. By interpreting the poles of the amplitudes in || as boundary 
bound states, we find a spectrum of masses that look like ( |2.10| ) but with a slightly differ¬ 
ent range of parameters. Considering then scattering of breathers off a boundary with a 
bound state does not give rise to any new poles beside (^] 2 |) and (| 2 . 9 |) , with in the latter 
case an extended range of values of N (for simplicity we do not give the relevant boundary 
S-matrices here). Therefore the complete boundary bootstrap is closed in principle. 


2.2. Integral representations of various S-matrices. 


For comparison with results obtained from regularizations of the sine-Gordon model 
it is useful to write integral representations of the boundary .S-matrices ( 12.ID , ( fT6j ) and 
( |2.7|) using the well-known formula 


iogr(z) 



dx _ T 
— e x 
x 


z-l + 


e -(z-l)x _ | • 

1 — e~ x 


Rez > 0. 


(2.13) 
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Suppose first that l<2£/7r<A+l and denote 


n* 


7T 


(2.14) 


where the square brackets mean the integer part of the number. For such values of £ 
there are n* + 1 poles ( j2.2| ) in the physical strip, i.e. the spectrum of excitations contains 
boundary bound states. Correspondingly, there is a finite number of T-functions in ( |2 . 1|) , 
(pip, ( [2.7| ) whose arguments have negative real part so that formula ( |2.13| ) is not applicable. 
Treating such T-functions separately, we obtain the following results: 


• d l 

- ! d« log 


p + (e) 

[R 0 (6) 


2A f + °° , (2X9x\ 

— — / ax cos - 

^ J-OO V 7T J 


X 


sinh(2£/7r — 2 n* — 2)x sinh(A — 2£ t /ir)x 


sinli x 


2 sinli x cosh Xx 


(2.15) 


■ d i 

-'Te log 


Ro(0) 


2A f +oc , f 2X9x\ 

— — / ax cos - 

oo V 7T J 


7T 


X 


sinh(A — 2£/ n)x — 2 cosh x sinh( A + 1 + 2n — 2£/n)x 


■ d ^ 

- l Te log 


Ro(0) 


2A /'+ 00 , / 2X9x\ 

— / ax cos - 

* J-oo \ 7T / 


2 sinli x cosh Ax 

sinh(2n* + 2 — 2£/7 t)x 
sinli x 


sinli(A — 2£/7r)x — 2 cosh x sinh(A + 1 + 2n — 2£/7r)x 
2 sinli x cosh Ax 


(2.16) 

(2.17) 


In the derivation of analogous representation for P there are no subtleties because the 
’’dangerous” T-functions cancel. We get 


■ d l 

- l Te log 


p-(0) 

L^o(0)J 


2A f + (2X9x\ smh(A — 2£/7t)x 

— / ax cos - ---——. 

7T oq \ it J 2smhxcoshAx 


(2.18) 


In the region 0 < 2^/ir < 1 where there are no poles and no boundary bound states in the 
spectrum, formula ( |2.18|) is valid, too. The expression for P + can be obtained from ( |2.15| ) 
by setting formally ub = n* + 1 = 0, which gives 


• d i 

■’Te log 


[Ro(0) 


2X f + °° A f2X9xX^ sinh(A + 2£/tt)x 

/ CLX COS I I ; ■ ■ ■ . 

7T J_oo \ 7T J 2 smh x cosh Xx 


(2.19) 
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Note that if 2^/ti > 1, the integral in ( 2.19 ) diverges. Finally, we complete this list by the 
following two expressions: 


■ d i 

- ‘Te log 


Ro(0) 


■ d i 

= ~ l Te log 




+ 


2A 


,,+00 , f2\6x\ 

H-/ ax cos 




R 0 (9) 
sinh(— — 2 n* — 2 )x 


sinli x 


For the integral representation of Rq see [BJ. 


2 cosh x sinh( ^ + 2N — X — l)x 
2 sinh x cosh Xx 

( 2 . 20 ) 


3. Exact solution of the regularized boundary sine-Gordon model. 

3.1. The XXZ chain with boundary magnetic field. 

The XXZ model in a boundary magnetic field 


n = 


7T — 7 
27rsiii7 


L — l 


Y («+1 + <Ti<Ti+i + 1 - !)) + h{°\ - !) + h'{a z L - 1 ) 


2=1 


, (3.1) 


was discussed in [| 10 |, where its eigenstates were constructed using the Bethe ansatz. As 
usual, these eigenstates TL\n) — E\n) are linear combinations of the states with n down 
spins, located at x±, ...,x n on the chain: 


I n) = ^/ (ri) (xi,...,x n )|xi,...,x n ). 

Consider for simplicity the case n — 1. The wave-function f^\x) reads [10 
/W(x) = [e~ ik + {h! - - (k -»• -k) = 


sinh |(zy + a) 
sinh |(ry — a) 


L—x 


sin 7 sinh ^(a + i'yH') 

— — (cr —> -a), 


sinh tt (z 7 — a) sin 4(7 + 7 -H 7 ) 


where we defined the new variables as in |nj: A = — cos 7 , k = /(ck, 7 ), 

h — i sin 7 


7 # = /(ry, — i ln(h — A)) = —7 — z In 


h + i sin 7 


(and similarly for H'), and 


/(a, 6 ) = —i In 


sinh k{ib — a) 
sinh |(i& + a) 


(3.2) 


(3.3) 


(3.4) 
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When h varies from 0 to +oo, 7 H increases monotonically from — tt — 7 to —7 according 


to ( |3.3|) if we take the main branch of the logarithm. 

Denote h t h = 1 — cos 7. This “threshold” value of h corresponds to 7 H = — 7 r; its 
meaning will become clear below. When h varies from — 00 to 0, 7 H increases monotoni¬ 
cally from —7 to 7T — 7. For the purposes of the present work we confine our attention to 
the region h, h! > 0 and choose 7 G (0, ^). Other regions in the parameter space can be 
obtained using the discrete symmetries of the Hamiltonian ( |3.1| ): a z —> —a z on each site 
or on the odd sites only. The parameter k in (|3.2|) is not arbitrary, but satisfies the Bethe 


equation 10 


j(2L-2)k (e ik + h- A)(e ik + h'- A) = 

(e~ ik + h - A)(e~ ik + h'- A) ’ 


or 


sinli | (a 


*7) 


2 L 


sinh |(a + iy) 


sinli — i^H) sinh i(a — zyiT) 
sinh + i^H) sinh h(a + i^H') 


= 1 . 


( 3 . 5 ) 


( 3 . 6 ) 


Note that the wave-function ( |3.2| ) depends on H implicitely through the solution of the 
Bethe equation ( |3.6| ) a(H,H'). Besides, one can multiply the amplitude (fT2j) by any 
overall scalar factor depending on a, L. H and H'. The Bethe equations in the sector of 
arbitrary n > 1 can be found in [ 110 . 


3.2. The Bethe equations for the inhomogeneous XXZ chain. 

The real object of interest for us is actually the inhomogeneous six-vertex model with 
boundary magnetic field on an open strip. The inhomogeneous six-vertex model is obtained 
by giving an alternating imaginary part At A to the spectral parameter on alternating 

It was argued in [; 6 j, generalizing known results 


vertices of the six-vertex model [ 11,12 


for the periodic case | 12 J that this model on an open strip provides in the scaling limit 
A, L —> 00, lattice spacing ->0a lattice regularization of ( |1.1|) , with (3 2 = 87 and a value 
of <po at the boundary related to the magnetic field. The reader can find more details on 
the model in the references; it is actually closely related to the XXZ chain we discussed 
above. In particular, the wave function can be expressed in terms of the roots ay of the 
Bethe equations |l(],|Il] : 


sinh | (a 3 + A — iy) sinh |(aj — A — zy) 
sinh 7j(otj + A + iy) sinh \ {otj — A + ly) 


sinh |(aj — i^H) 
sinh |(aj + i^H) 


x 


sinh ^ (aj 


i^H') 


sinh + i'yH') 


n sinh i(a 


a m - 2iy) sinh t ( a? + a m - 2iy) 


m.^j 


. sinh \ {cxj — CL m + 2iy) sinh + dm + 2iy) 


( 3 . 7 ) 



































By construction of the Betlie-ansatz wave function, a 3 > 0. Note that the solutions of 
([L7|) aj — 0, in should be excluded because the wave function vanishes identically in this 
case. The analysis of solutions of the ( |3.7|) is very similar to the case of the XXZ chain 
in a boundary magnetic field. We consider the regime 0 < 7 < 7 r/ 2 , which falls into the 
attractive regime 0 < f3 2 < An in the sine-Gordon model ( | 1 . 1 | ). We set hereafter 7 — n/t 
and for technical simplicity restrict t to be positive integer. In the limit L —> 00 this 
constraint implies that in the bulk only the strings of length from 1 to t — 1 are allowed, 
together with the anti-strings. 

Taking the logarithm of eq. (T7), one obtains: 


L [f (aj + A, 7) + f( a j ~ A, 7)] + 7 H) + f(oij , 7 H') 

= 2 t dj + Y [f( a j ~ «m, 27) + f(aj + a m , 27)] , 

m^j 


(3.8) 


where lj is an integer. We also recall the formula for the eigenenergy associated with the 
roots a 3 10 

9( 7T — - 

(3.9) 


E = ——— + A . t) + /'(«J - A, 7 )]. 


n 


3.3. Thirring model with boundary. 


Like the bulk sine-Gordon model is a bosonized version of the bulk massive Thirring 
model, one can expect that the boundary sine-Gordon model is a bosonized version of 
the Thirring model with certain boundary conditions. The quickest way to identify this 
boundary Thirring model is to use the Bethe ansatz equations ( |3.7|) . Write the most 
general t/(l)-invariant boundary interaction 


H t = dx[-iipf ipi x + i^ V’ 2 z + m 0 ^ + ip 2 + m 0 ^i + 2#oV , i~ 

Jo 

+ Y ai ^i ^ (°) + a '^t^3 (A) ■ 


(3.10) 


13 


13 


The entries of the 2x2 matrices A — A' = {ah} can be determined up to one arbi¬ 

trary parameter ip by the hermicity of Hr and the consistency of the boundary conditions 
(detA = 0). For the left boundary, the matrix A looks like 


1 /e -< * 1 \ 

2 sin <j> \ 1 J 


(3.11) 
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and the boundary condition reads Vh(0) = — e*^ 2 ( 0 ) (similarly for the right boundary). 

To find the eigenvectors of the Hamiltonian ( |3.10| ), Hr T = E\ H, one can use the same 
wave-functions as for the bulk Thirring model (TT|| . and modify them by analogy with the 
example of XXZ chain in a boundary magnetic field [10]. This way one gets the equations 
for the set of rapidities a 3 : 

cosh + i(f>) cosh 7 ? (ay + 


0 2 imoL sinh a.j 


cosh ^(aj — icj)) cosh ^(aj — icj)') 

-|-j- sinh |(aj — a m — 2iy) sinh + a m — 2 i'y) 
sinh \ (a : j — a m + 2 i'y) sinh k( a j + a m + 2 i'y) ’ 


(3.12) 


x 


m^j 


where 7 is related to go in the usual way B These equations look quite similar to 
( 13.71 ). Tim mapping can be made complete by taking in ( [3.71) the limit A —> 00 with the 
identification mo = 4 e _A siii 7 . 

The derivation of these equations is rather cumbersome, therefore to illustrate the 
procedure we comment on the simplest case of one-particle sector, which is nevertheless 
sufficient to obtain the form of the boundary terms in ( |3.12| ). We make an ansatz T = 
Jcf dyx^iy)^ (y)|0), where A is the spinor index, x{y) is the wave-function and |0) is the 
unphysical vacuum annihilated by 'W\. 

The equation HyT = E\\> reduces to (a^ are the Pauli matrices): 


d 

-At 3 — x + m 0 a ix + Ax6(x) + A'x5(x — L) — Ex. 


We look for the solution of ( |3.13| ) in the form 

,-<*/2 


Ai 

X2 


= a(a) 


W 2 


0 im qx sinh a. 


a 


-a) 


e “/ 2 

3 —a/2 


j—im qx sinh a. 


(3.13) 


(3.14) 


Substituting it into ( |3.13| ) we get E — m 0 cosh a and, besides, two boundary conditions to 
be solved. The first one, at x — 0, determines the form of the factor a(a) — cosh A (a — i(j)), 
while the second one at x — L gives rise to the Bethe equation 


2zmo-Lsinha _ 


cosh i(a + i(j)) cosh h(a + icj)') 


cosh | (a — i(f>) cosh |(a — icj)') ’ 

which determines a. Comparing the Bethe equation ( |3.7|) with ( [ 1 . 1 2 | ) and using the relation 
between £ and El obtained below in section 5 we find the relation between the boundary 
parameters (j) and tpo in the Hamiltonians ( |3 . 1 (( ) and ( | 1 . 1 | ) respectively: 

cj) = (3<p 0 - /5 2 /8. 
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Thus, the integrable boundary condition for the f/(l)-in variant boundary Thirring model 
reads: 


^ 2 (0) = -e l/?2/8 -^°^i(0). 


(3.15) 


It would be interesting to obtain the result (|3.15|) directly from the Hamiltonian ( |1 . 1| ) 
using an extension of the Coleman-Mandelstam bosonization technique to the case with 
boundary. However, to our knowledge such an extension has not been developed yet. The 
naive application of the known Coleman-Mandelstam “bulk” formulas doesn’t give the 
factor e l P / 8 in (p,15|) , which seems to be some kind of “boundary anomaly”. 


4. Solutions of the Bethe ansatz equations with boundary terms. 

As is well known in the case of the bulk Thirring model or equivalently the peri¬ 
odic XXZ chain, the bound states are associated with various types of solutions of the 
Bethe ansatz equations involving in general complex roots |T3[. By analogy, we expect 
the boundary bound states to correspond to new solutions made possible by the boundary 
terms. 

Consider first the example of the XXZ chain as given in section 3.1. Since our goal 
is to study purely boundary effects, we will look for the solutions of the Bethe equations 
that give rise to a wave-function localized at x — 0 or x — L and exponentially decreasing 
away from the boundary. 

The states described by such wave-functions will be referred to as the “boundary 
bound states” below. For this, one should have a purely imaginary in (|372]). We consider 
here the limit of L large, when the left and the right boundaries can be treated inde¬ 
pendently and the overlap of the corresponding wave-functions is negligibly small (for the 
physical applications it is necessary to take the scaling limit anyway). In the limit L —> oo 
it is easy to check that there are two such solutions to ( |3.6|) : a — ia o = H, H') 

and a — ia' 0 — —i^H' + ie'(L, H, H'), where £ ~ exp(—2 kL) and we defined k > 0 as 


e — 


sm 


[~lH - 7 ) 


sin ±(—7# -(- 7) 


(similar relations are assumed for e', k'). Solution a' 0 gives a wave-function (|3.2|) localized 
at x — L: /W(x) ~ e~ K ( L ~ x \ Solution q : 0 gives a wave-function localized at x — 0, 
~ e -Ka: , provided we renormalize the wave-function (372): 

1/2 


/ (1) - / (1) 


sinh -(a — i'yH) sinh - (a + i'yH) 

2 2 


(4.1) 
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In the special case H — H' there is only one proper solution a — iao — —i^H + is(L, H ) 
with £ ~ exp( — kL). The wave-function ( |3.2|) behaves as the superposition of the ’’left” and 
the ”right” boundary bound states, ~ ( e~ KX + g-'h 1 '-*)). Note that the boundary 
bound state appears in the above example only when the boundary magnetic field is 
large enough: namely, h > hth ! • This follows from the fact that a should be such that 
0 < a 0 < n. 

Now, consider the equations for the inhomogeneous model ( |3.7| ). The basic boundary 
1-string solution to ( |3.7| ) is a — iao — —i'yH + ie, provided that 0 < ao < 7r. This solution 
is possible due to an argument very similar to the one used in the bulk: as L —> oo, the 
two first terms of ( |3.7| ) decrease exponentially fast, while the third increases exponentially 
fast, and e ~ exp(—2 kL) with 


sinh 2 4 + sin 


2 ap -7 


e — 


sinh 2 4 + sin 2 


(4,2) 


Recall that for the bulk problem when there is no boundary term, the right hand side of 
( [3.71 ) would have to decrease exponentially, forcing the existence of a “partner” root at 
a — 2ry. 

One can construct similarly boundary n-strings which consist of the points iao, ' l(y o + 
2zy, o + 2 i{n — l)y (see figure 1). By convention n — 0 means there is no boundary 
string, that is all complex solutions are in the usual bulk strings. The possible values 
of n are restricted by the fact that the upper point of the complex should be below in: 
max(n) = | " .4 >0 ] + 1, where the square bracket denotes the integer part. To show that 
the boundary n-string is indeed a solution to (|3]7]), we introduce infinitesimal corrections 
£j to the positions of the points of complex |TH|. Taking the modulus of both sides of 


2 More generally, the criterion of existence of boundary bound state solutions allows us to 
determine threshold fields for any A. For this, let us examine (|3.5| ). The parameter k is defined 
modulo 27 r, therefore we restrict it to lie withing k G ( 0 , 27 r). Two possibilities k = ia and 
k = 7r + ia, where a > 0, lead to two different threshold fields, determined by the fact that the 


denominator in (T5) should vanish: 


hth — A + 1, 


= A - 1, 


and the regions with boundary bound states are h < A — 1, h>A+l. When A > 1, there 
are two different threshold fields, in agreement with the results of Jimbo et al |14[|. In the region 


| A | < 1, discussed here, there is only one threshold field h 


(i) 
th ’ 
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( p~7| ) with aj — ia o + 2ik^ and multiplying equations for k = 0, ...,n — 1 we obtain 
exp{—2L(«i + K 2 + ... + K n )} ~ £1, where E\ denotes the correction to the point iao . The 
behavior of the remaining £k follows from £ i by recursion. For example, for the 2-string 
£2 is given by \e\ — £2] ~ exp(—2 Lk 2 ) 

Additional boundary strings can be obtained by adding the roots ia s below iao so that 
ia s — iao ~ 2is7, with s = 1,2, ..., N (see figure 2). Together with the boundary n-string 
above ao , they form the complex which we call boundary (n, N)- string. To analyze the 
existence of such complexes as the solutions of ( |3.7| ) we introduce as before the infinitesimal 
corrections £ s to the roots a s , where now s — n,n — 1,1, —1, —2,— N. Then, the 


equations ( |3.7| ) with ay = ia s tell us that the range of N should be 


«0 jy- < 7T + cto 


27 


27 


(4.3) 


In other words, the inequality ( | 4 . 3 |) states that the lowest root of the boundary string 
should be below the axis Ima = 0 and above the axis Ima = —m. Another constraint 
follows if we multiply the equations ( | 3 . 7 | ) for all the roots of boundary (n, IV)-string. This 
gives exp(—2 LJ2 K s) — £ i- So, one should have k s > 0. The latter sum can be 
easily evaluated if one uses the expression ( (472|) simplified in the limit A —> 00: n — 
4e -A sin 7 sin ao- The constraint obtained in such a way forces the number of roots above 
Ima = 0 axis in the boundary string to be greater than the number of roots below Ima = 0. 

We have not been able to find any reasonable additional solution to the Bethe ansatz 
equations. The two sets of boundary strings we have encountered appear to be in one 
to one correspondence with the boundary bound states identified in section 2 using the 
bootstrap approach. To clarify this identification we now compute related masses and 
S-matrices. 


3 Note that associated with each boundary n-string there is also the solution to (p.7|) obtained 
by complex conjugation of all c/s. The existence of such a ’’mirror image” is the consequence of 


the symmetry of equations (3.7) and it is of no importance to physics. In the bulk case, it is easy 


to show [ 16] that all solutions are invariant under complex conjugation, but this result does not 
hold here. In fact, a solution which both the boundary n-string and its mirror image would lead 
to a vanishing wave-function. 
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5. S matrices and bound state properties from the exact solution. 

5.1. Bare and physical Bethe ansatz equations. 

The ’’bare” Bethe equations follow from taking the derivative of (|3.8| ). Defining 
2L(pk + Pk)dct to be the number of roots in the interval da, one obtains coupled inte¬ 
gral equations for the densities of strings p\, ...,pt -i and anti-strings p a : 


2n(p k + Pk) = \v'h ~ fL * Pa ~ fki * Pl+ 2 T^ Uk ~ ~ U '^hk) 

i=i 

2n(p a + Pa) = ~\p'a + /' * Pa + ^ fal * Pi + + u'f'^) 


t -1 


(5.1) 


1 = 1 


where * denotes convolution: 


f*g(a)= d(3f(a — f3)g(/3). 

J — OO 

These densities are originally defined for a > 0, but the equations allow us to define 
Pk(—a) = pk(o:) in order to rewrite the integrals to go from — oo to oc. If we totally neglect 
the boundary terms (terms ~ L~ l ) in ( |5 . 1|) , we will end up with the same equations as 


for the periodic inhomogeneous six-vertex model [11]. The various kernels and sources in 
flO|) are defined as follows: 


Pa{ot) = f(in + a + A, 7) + f(in + a - A, 7), 

Pk(a) = f(a% + A, 7) + f(ai - A, 7), 

OLi 

where the sum in the last expression is taken over the rapidities of the bulk k-string root 
centered on a. 

The kernels fki are the phase shifts of bulk k-string on bulk /-string obtained by 
summing (|3.4|) over the rapidities of string roots. The boundary terms are: 


= -2/'(2a, 2y) -/'(« + in, 7#) - f\a + in, 7 H') - 2n5(a), 
u k{ot) = J^[ 2 /'( 2 ai, 27) + f'{a.i, 7 H) + f'fa, 7 H')\ - 2n5(a), 

Oii 

(the sum above is over the roots of bulk k-string centered on a), 

fn L,R \ a ) = '^2f(in + a-a i ,2'y) + f(in + a + a i ,2'y), 
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and Qij denotes the rapidities of the roots in the boundary n-string. 


ffy R \ a ) = 2 to+ f( a j + 2 ^)’ 


where cq denotes the roots in the boundary n-string, while ay denotes the roots in the 
bulk k-string centered on a. The parameters u>,u/ are equal to 1 or 0, depending on 
whether the boundary string is present or not . In our case, 0 < 7 < 7 t/2 , the ground 
state of the periodic inhomogeneous XXZ chain is filled with anti-strings. The physical 
Bethe equations are obtained [|17|,[|T8| by eliminating the ”non-physical” density p a from 
the right-hand side of ( p.l| ). This is done simply by solving for p a in the last equation in 
( pTTp and substituting it into the others. The result is 

2 ^(pk + 14 ) = Ij4 + \ 2 I-f * p « + 


+ A * - X (fL + ^4) * Pi + 


1=1 


2 TT-f' 

t -1 


2 n-r 

1 2ir 

2 2 V-f 2n — f' ~ " ra 1 2n — f 

u a + + ui'f'S R) 


(5.2) 


9 / , .fcx _ 1 27T P'g _ f . 9 _ n h , ST Pal * 9 ™ +2 

^ Pa Pa ^ 2 2 t r - P 2 tt - f' Pa f-' 2 tt - f' Pl + 2L 

where 


u, 


n,n ' ;a? 


Ur, 


= 2tt- 


27T - /' 

Un,n'\k =U k - ufffi - U/- f' ak * U n , n ’- a /2Tl, 

and different products (ratios) of kernels are defined through their Fourier transforms. 


(5.3) 

(5.4) 


5.1?. The mass spectrum of boundary bound states. 

We assume at first that the ground state is built by filling up the Dirac sea with 
anti-strings, as in the case of the periodic XXZ chain. We will see below that this is not 
always true. The presence of the boundary strings in the Bethe equations deforms the 
distribution of roots and modifies the density of the Dirac sea p a by a term 5p a /2L of 
order L -1 . With the boundary n-string, the Bethe equation for the density of the Dirac 
sea particles p a is 

1 1 f + °° , 1 

2 Pa( a ) - 2 l u a( a ) = J f (<x-P)pa(P)dP-2np a (a) + —f n (a), (5.5) 
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where f n was defined above. Subtracting from (5.5) the equation for the density of the 
Dirac sea alone, 


2 Pa( a ) ~ T^M") = 


- 4-00 


/'(a - (3)p a {(3)df3 - 27rp a (a), 


one obtains the equation for 5p a : 

f+OO 


0 = - 


f'(a - p)5p a (P)dP + 2tt 5p a (a) - f^a) 
5pa = 2 L(p a - p a ). 


(5.6) 


(5.7) 


The solution to (5/7) can be written in terms of the Fourier transform 5p a {k) — 
f dae lka Sp a (a) as follows: 

fnW 


dpa(k) = 

2 tt - f'(k) 

For the boundary n-string iao + 2zy s, s — 0, 1,..., n — 1 we obtain 

4 coshy/c sinh ny/c cosh(o:o + yn — y)/c 
sinli 7T k 

2 cosh y k sinh ny/c cosh(ao + yn — y)fc 


fh( k ) = ~ 27r 
= 


where we used 


/'(/c) = 2 tt 


sinhy/c cosh(7r — y)fe 
sinh(7T — 2y)/c 


(5.8) 

(5.9) 
(5.10) 


sinh 7r k 

Expressions (|5.9| ), ( |5.10| ) are valid for the n-strings with n = 1, 2,..., • For the longest 

n-string with n = +1 = n* + l the Fourier transforms f' n . Sp a differ from ( |5.9|) , ( |5. 10Q : 

2 sinh(7r — 2y)A: cosh(cto + 2yn* — n)k 


K( k ) =2n- 


- 2tt 


sinh irk 

4 cosh y k sinh n*y/c cosh(cro + yn* — y )k 


(5.11) 


bPa(k) 


sinh irk 

sinh(7T — 2y)fc cosh(«o + 2y n* — n)k 
sinhy/c cosh(7r — y )k 

2 coshy/c sinhn*y/c cosh(cK 0 + yn* — y )k 


(5.12) 


sinh y/c cosh(7T — y )k 

The conserved U( 1) charge in the boundary XXZ chain is the total projection of the spin 
on the z-axis. In the thermodynamic limit the charge of the boundary n-string with respect 
to the vacuum is determined by [|T0j : 


0 + OO 


o + OO 


Qn n T 


2 Lp a da — 


2Lp a da — n + 


1 


o + OO 


Sp a da = n + -5p a (0). (5.13) 
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Using ( |5.10 ), we obtain for the n-string Q n = 0, and for the longest boundary string Eq. 
( |5.12|) yields Q* = ix / 2'-/. Similarly, the mass of the boundary strings in the thermodynamic 
limit according to (fU9[) is given by 


-+oo 


r> + 00 


" + OO 


m n = h n + / 2 Lp a h a da- / 2 Lp a h a da = h n + - / h a 5p a da, (5.14) 


where the expression for h a is 


- 2(tx — 7 ) . , . . 2 sinh 7 fccos Ak 

ha(k) = -- -Pa = —4(tt - 7)- 


7T 


sinli rxk 


and the soliton mass 11 


TO = 4e 2(n -y) _ 


We obtain in the limit A —> 00 


TOrj = to 


sin (2n — 1 — H) + sin (H + 1) 
zA 2 A 


7T 


to* = to sin — (iJ + 1 ). 
zA 


(5.15) 


(5.16) 


Since the parameter H varies in the interval —A — 1 < H < -1, the mass of the longest 
string to* ( 5. 16|) is always negative, while the other boundary strings have positive masses 
( |5.15|) . This means that the vacuum we have been working with is an unstable one in the 
region — t < H < — 1 (h > hth)- To cure the situation we define a new correct ground state 
by attributing the longest boundary string to the Dirac sea. The boundary excitations are 
obtained by succeessive removing of particles from the top of the longest boundary string. 
The charge and mass of such excitations with respect to the correct ground state are given 
by 

7r 7r 

Qn = - —, TO n = TOCOS — (A + 1 + H - 2n), n — 0,1,..., n*. (5.17) 

Z 7 zA 


Note that the number of excitations ( |5.17| ) is equal to the number of particles in the longest 
boundary string, n* + 1. The charge of such boundary excitations is equal to the charge of 
the hole in the Dirac sea. We identify a hole with a sine-Gordon soliton, and the boundary 
excitations described above, with the boundary bound states \/3 n ) ( |2.2| ). Their masses and 
charge ( |5.17|) and the counting coincide provided that 


t + H + 1 = 

TC 


(5.18) 
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and the lattice charge Q is properly normalized. This expression is in fact valid for all 
values of h > 0. The authors of 0, deriving this relation in the region h < hth, obtained 
a different expression because they used different branch of logarithm in ( |3.3| ) . 

In the above discussion we considered the boundary bound states related to one of 
the boundaries (say, the left one). In principle, one should include into the ground state 


the longest boundary string ia' 0 + 277 /, / = 0, 1,..., 


t+H 1 

2 


, corresponding to the right 


boundary as well. The energy of the excitations due to both boundary strings is a super¬ 


position of energies of the form (|5.17|) . When H — II'. these two boundary strings overlap 
and the usual Bethe wave-function vanishes. However on physical grounds we do not ex¬ 
pect anything special to happen when the boundaries are identical. So, in such a case one 
should use as a wave function a properly renormalized version of the limit H —> H' of the 
usual Bethe wave function. 

When the magnetic field varies, the above picture indicates a qualitative change in 
the structure of the ground state at values H = —/,. —t + 2,—t + 4,.... At these values, the 
mass of the bound state with the highest mass approaches the soliton mass and it becomes 
unstable. As discussed in [Q for the Ising case, this decay corresponds to large boundary 
fluctuations that propagate deeply into the bulk. 

The mass of the boundary (n. TV)-si ring with respect to the correct vacuum can be 
calculated analogously. The result is: 


mn,N — m cos( 


—) + mcos(— 


2n+l , , 2N- 

__ x) _ m c OS (- + — 


*■)> 


(5.19) 


where we used ( |5.18|) to express H interms of £. This result is rather confusing to us, 
because the above mass does not correspond in general to one of the bound state masses 
found in the bootstrap apporach. It can be considered as a sum of such masses, hinting 
that the (n, N ) string describes actually coexisting bound states, but the corresponding 
boundary S-matrix does not allow such an interpretation. We are forced (but see the 
conclusion) to consider that only the (n, 7V)-strings with n — n* + 1 occur, that is the 
physical excitations are built by adding roots to the ground state configuration below iao . 
The charge and energy of such excitations with respect to the correct vacuum is given by 


Qn — 


7r 


7r 


2y 2y 


= 0, rriN — m cos( 


7r 


£ _ 

A 2A 


£ 2N - 1 

) — m c °s("^ H ^ tt) 


(5.20) 


These coincide with the charge and mass of the boundary bound state |5 n= o,Ar > (Q). 
The range of N (|4.3|) agrees with the range of corresponding parameter in ( |2.9| ). 
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5.3. Boundary S matrices. 


It remains to check that the boundary S-matrices obtained above by the bootstrap 
approach coincide with those of the lattice model. To extract the boundary S-matrices 
from the Bethe equations we will follow the discussion of ||. Briefly, the idea of the method 
is the following. The physical excitations of the lattice model in the limit A —> oo can be 
thought of as relativistic quasi-particles with rapidities 6i. The iutegrability implies that 
the set {9i} is conserved. Moreover, if the scattering matrices are diagonal, each particle 
preserves its rapidity. Assuming that this is the case, the quantization of a gas of A f quasi¬ 
particles on an interval of lenglit L results in the integral equations for the set of allowed 
rapidities ||: 

p 1 

2ir(p b + p^) = m b cosh d + ^ (p bc * p c + — 0 & , (5-21) 

C= 1 

where subscript stands for the type of particle, and 


<Pbc(0) = -i^ln S bc (6) 
Qb(e) = -i-^\nR b(L) 


p ( 9 ) - i ^ \nR b 0 ( , R) {6) + i— In S bb {29) - 27tS(9). 


d 


(5.22) 


Equations ( 5.2 1| ) should be compared with the physical BE ( |5 . 2| ) , which gives bulk and 
boundary S-matrices. We will confine our attention to the boundary S-matrices only, 
keeping track of those terms in ( |5.3|) , ( [All) , ( |5.22|) , which depend on the boundary magnetic 
field (the field-independent terms contribute to Rq and their agreement has been shown in 
[(|). The discussion for the left boundary is completely parallel to that of the right one. 
Also, it is sufficient to consider only b =solit.on and b =anti-soliton in ( |5.21|) . We identify a 
hole in the anti-string distribution in ( |5.2|) with a soliton in ( 5.21 ), and (t — l)-string with 
an anti-soliton. Below we give explicit expressions only for the kernels in ( [5 . 2|) which we 
need for our analysis. The other expressions are listed in |J. Suppose first that h < h cr 
(—t — 1 < H < —t). This corresponds to the case without boundary excitations in the 
spectrum, £ < 7 t/ 2. Choose u> = u/ = 0 in (|5.2|) . Then 


ww = 2 W lh(2 T7 g W 

sinh7rfc 


(we omitted the I^-independent terms and iT-clependent ones), 


fj{L) = 2?r siiih(27r + 7 g)fc 

2 sinhy/c cosh(7r — y)/e 
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Using ( 5.22 ) we compare this expression with ( 2.19|) (recall that the rapidity a should be 
renormalized a —> 0 — ta/ 2A) and find complete agreement under the identification ( |5.18|) . 
Similarly, one can use 


u 


(L) _ 


t -1 


—2tc 


sinh(7r + 7 H)k sinh(7r — 7 )k 
sinh 7 rk sinli 7 k 


+ . .., 


ft(L) = 2rr sinh ( 2 + H h k 

1 2 sinh 7 A; cosh(7r — 7 )k 


+ ... 


to compare Ut -1 with ( |2.18|) and obtain agreement as well. Next, suppose that h > hth 
(—t < H < —1). To obtain the boundary S-matrices for scattering on the ground state 
|0)s set u — u' — 1 and choose the boundary string to be the longest string, n — n* + 1 
in (U). Then, using ( pT)) , = -/* and 


u 


(L) = 27r sinh 1 Hk 


sinh 7 rk 


u 


( -l = - 27T 


+ . . . , 


sinh(i? + 2[^-]) 7 fc 
sinh 7 k 


+ ... 


(5.23) 


we obtain 


U, 


(L) 

n* + l;a 

2tt 


sinh 7 Hk 


2 sinh 7/c cosh(7r — 7) k 


+ 


sinh(7r — 27)^’cosh(iJ + t — 2n*)^k 
sinh 7/c cosh (77 — 7) k 


2 COSI17A’ sinh n^k cosh(if — n* + 1)7^: 


sinh 7/c cosh (77 — 7 )k 


+ ..., 


U. 


(L) 


n* +1 


- rt L ) 




„ sinh(/7 + 2 [i=ii]) 7 * : 

— 271 -—-;-b • • • , 


n * +1;a sinh 7/c 

which agrees with ( |2.15|) , (12.18 ) under the identification ( 5.18 ). Note that the last relation, 


which follows directly from (|5.3|), (|5.4|) and ( |5.23|) , is valid also for U n - a and U n -t- 1 with 
any n. In the same manner one can calculate the boundary S-matrices for scattering 
on the boundary n-strings and check that they indeed coincide with ( |2 .1 fi| ) , ( |2.17| ) under 
the condition (|5.18|) . For this, one needs to take u = u' = 1 in (|5.2|) and use ( |5.9f ), 
fn-t-i — ~fn ■ Finally one can compute also the boundary S-matrix for the scattering on 
the (n* + 1, IV)-strings, again in agreement with the bootstrap results. 
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6. Conclusion 


The question of boundary bound states even in the simple Dirichlet case appears rather 
frustrating: using the XXZ lattice regularization or equivalently the Thirring model, we 
have only been able to recover the f3 n and h n =o,iv boundary bound states. A way out is to 
consider solutions of the Betlie ansatz equations made of an (n, N ) string superposed with 
the n* + 1 string that describes the ground state. This is not, allowed in principle in the 
model we consider because the Bethe wave function vanishes when two roots are equal. 
However, putting formally such a solution in the equations gives the masses of the 5 n ,n 
states and the S-matrix also agrees with the bootstrap results! But the meaning of this is 
not clear to us. 

Finally let us mention that one can calculate the ground state energy in the ther¬ 
modynamic limit by solving the equation (|5]5|) for the ground state density and using 


(0): 


r»+oo 


E gr 


2 Lp a (a)h(a)da. 


As a result we get the combination E , 


gr 


sine-Gordon ground state energy [19 


E bu ik + E boundary , where E buik is the well-known 




Lmr 


tan 


7T7 


2(tt - 7) 


and E boun dary is the contribution of the boundary terms: 


Ei 


m 


boundary\ 


• (J/+2) 7 7t 

2( +1 + cot * 


sin 


2(tt-7) 


4(7T — 7) 


We see that the ground state energy of the boundary sine-Gordon model is a smooth 
function of the boundary magnetic field for the whole range of h in the XXZ regularization, 
hence of <po- The changes in ground state structure do not affect E. as is expected since 
in such a unitary model there is no (one dimensional) boundary transition. 

The finite size corrections to the ground state energy themshelves (the genuine Casimir 
effect) can be computed using the technique developed in [|19[. It is also interesting to 
consider the inhomogeneous 6-vertex model with an imaginary boundary magnetic field 
insuring commutation with U q sl( 2) [[20J . This should presumably lead to a solution of 
minimal models with integrable boundary conditions. We will report on these questions 


soon. 
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Figure Captions 


Figure 1: the first type of boundary string. In the ground state, the boundary string of 
maximum allowed length is occupied. 

Figure 2: the second type of boundary string. 
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